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ABSTRACT 

Investigating the dynamical evolution of dust grains in proto-planetary disks is a key issue to 
understand how planets should form. We identify under which conditions dust settling can be 
constrained by high angular resolution observations at mm wavelengths, and which observa- 
tional strategies are suited for such studies. Exploring a large range of models, we generate 
synthetic images of disks with different degrees of dust settling, and simulate high angular res- 
olution (~ 0.05-0.3") ALMA observations of these synthetic disks. The resulting data sets are 
then analyzed blindly with homogeneous disk models (where dust and gas are totally mixed) 
and the derived disk parameters are used as tracers of the settling factor. Our dust disks are 
partially resolved by ALMA and present some specific behaviors on radial and mainly vertical 
directions, which can be used to quantify the level of settling. We find out that an angular res- 
olution better than or equal to ~ 0. 1" (using 2.3 km baselines at 0.8mm) allows us to constrain 
the dust scale height and flaring index with sufficient precision to unambiguously distinguish 
between settled and non-settled disks, provided the inclination is close enough to edge-on 
(i ^ 75°). Ignoring dust settling and assuming hydrostatic equilibrium when analyzing such 
disks affects the derived dust temperature and the radial dependency of the dust emissivity 
index. The surface density distribution can also be severely biased at the highest inclinations. 
However, the derived dust properties remain largely unaffected if the disk scale height is fitted 
separately. ALMA has the potential to test some of the dust settling mechanisms, but for real 
disks, deviations from ideal geometry (warps, spiral waves) may provide an ultimate limit on 
the dust settling detection. 

Key words: Stars: formation — stars: circumstellar matter — ISM: dust 



1 INTRODUCTION 



Grain growth and dust settling are two key ingredients in the plane- 
tary system formation process. Recent observational evidences sug- 
gest that ISM dust grains start to grow in the early phase of star 
formation, as soon as dense pre- stellar cores begin to form. Theory 
and numerical simulations predict that in Class II proto-planetary 
disks, the dust orbiting the Pre-Main-Sequence (PMS) star contin- 
ues to grow but also very quickly settles along the mid-plane in 
typic a l characteristic time of a fe w 10 4 yrs foullemond & Dominikl 
12004 iFromang & Nelsonll2009h . The growth is the first step to- 
wards the formation of even larger solid bodies, which ultimately 
culminate with planetary embryos. Settling will speed up this pro- 
cess by favouring grain collisions, firstly by increasing the relative 
vertical velocities, as settlin g acts differently in fu nction of the dust 
dynamic properties (see e.g. lBirnstiel et alj2010l) . and secondly by 
concentrating dust close to the midplane. On top of that, a high 
dust to gas ratio in this area, can affect the gravitational stabil- 
ity and control the initia l step of the formation of planetesimals 
dGoldreich & Wardlll973h . 

Quantifying the dust evolution process is a complex prob- 



lem since the two physical processes (grain growth and dust set- 
tling) are simultaneously shaping the disk. The big grains are 
expected to fall down relatively quickly to the mid-plane while 
only small grains, r eflecting the stellar light (burrows et al.ll 19961 : 
lRoddieretaLlll996h . should remain located on the disk surface, at 
3-5 gas scale heights. 

At a radius of 100 AU from the central star, typical hydrostatic 
scale heights range between 10-20 AU or ~ 0.1" at the distance of 
the nearest low-mass star forming regions (D~ 140 pc). Therefore, 
observing settling requires both the most sensitive and the most 
resolving astronomical facilities. 

Some evidence of dust settling has been obtained from 
studies using Near-I nfrared (NIR) maps obtained by the HST 
fouchene et alj|2003h. o r by the analysis of th e Silicate band at 
lOum dPinte et all 120081 : ID' Alessio et all 1200 lb . IR observations 
only characterize grain growth for small particles with sizes a ~ 
0.1 — 10 /xm, as images at wavelength A are mostly sensitive to 
particles of size a ~ X/(2ty). Moreover, as the dust opacity in the 
NIR is still quite large, the particles we observe are necessarily lo- 
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cated high above the disk p lane, typically around 3-5 scale heights 
dChiang & Goldreichlll997h . 

Contrary to IR, the moderate opacity of the mm/submm do- 
main should probe material throughout the disk structure. The early 
bolometric observatio ns of envelopes and disks around young stars 
(iBeckwith et al.|[l99Qb indicated that both the dust absorption coef- 
ficient k u and its spectral index (3 at mm wavelengths have evolved 
compared to the ISM dust. However, only spatially resolved obser- 
vations could alleviate the ambiguity left by the possible contribu- 
tion of the inner optically thick core. Furthermore, contamination 
of the long wavelengths (longer than 4 mm) flux density by free- 
free emission can be substantial and should be removed fo r proper 
determination of the spectral index dRodmann et aljEoOol) . Using 
the VLA (at 7 mm), PdBI and OVRO to pr obe the dust prop erties 
and the dust disk surface density in CQ Tau jTesti et al.1 (120031) con- 
cluded that particle have grown up to sizes as large as ^ 1 cm. Sim- 
ilar results were obtained on larger samples in p Oph (with ATCA) 
and in Taurus-Auriga by iRicci et al . (2010b r a). The overall grain 
growth in proto-pla netary disks thus s eems a well establish fact. 

More recently, iGuilloteau et al.l (120111) performed a high an- 
gular resolution dual frequency study of disks in the Taurus -Auriga 
region with the IRAM array. Apart from disks with inner holes such 
as LkCal5 JPietu et alj2006h . all sources observed with sufficiently 
high angular resolution (0.4-0. 8 ") exhibit steeper brightness gradi- 
ent at 3 mm than at 1.3 mm. This is the signature of an evolu- 
tion of the dust spectral index with radius, with smaller f3 values 
near the central star. The inner part of disks, up to 60-80 AU, ap- 
pears dominated by large particles leading to a spectral index /3 
below 0.5 between A = 3 and 1.3 mm while beyond 100 AU, /3 
reaches value consistent with ISM-like grains (1.7). This consti- 
tutes the first observational evidence of radial variations in dust 
properties, and the characteristic transition radius between small 
and large grains is consistent wi th recent models of dust evolution 
in disks bv lBirnstiel et all bOlOh . 

In this paper, we go one step further and study the impact 
of dust settling on the disk imaging at mm wavelengths, in order 
to define adequate observational strategies to constrain this phe- 
nom enon with ALMA . For this purpose, we utilize the code DISK- 
FIT JPietu et alJl2007h . which has been upgraded to take into ac- 
count the dust se ttling. The ALMA simulator developed at IRAM 
(IPetv et al.ll2002h is then used to generate realistic ALMA datasets 
within the wavelength range 0.5 to 3 mm. Finally, we analyze these 
synthetic observations (pseudo-observations) as real data assum- 
ing a vertically uniform dust distribution in order to find out robust 
criteria of dust settling. We also explore some hidden degeneracies 
which may bias our estimate of the dust properties. We then discuss 
what would be an ideal ALMA observation. 

Our dust disk models are described in Section |2] Section 
[3] presents the ALMA predictions (pseudo-observations) and the 
method of analysis. We then discuss in Section 01 the implications 
of our results. 



2 MODEL DESCRIPTION 
2.1 Disk Model 

As in lGuilloteau et alj (120111) . we assume a simple parametric disk 
model. In Model 1, the gas surface density is a simple power law 
with a sharp inner and outer radius: 



S 9 (r) 



foriimt < r < .Rout- 

In Model 2, the density is tapered by an exponential edge: 



(2) 



S 9 (r) = So (-^-) " 7 exp (-(r/i? c ) 2 ^) . 

Note that Model 1 derives from Model 2 by simply setting R c 
oo and p = 7 in the above parametrization. Model 2 is a solution of 
the self- similar evolution of a viscous disk in which the viscosity 
is a power law of the radius (with constant exponent in time 7) 
(iLvnden-Bell & Prinddl 19741) . 

The kinetic temperature in the disk mid-plane is also assumed 
to be a power law of the radius: 



T fc (r,^ = 0) =T (r/Ro) 



(3) 



We assume that grains and gas are fully thermally coupled, so that 
the dust temperature Td us t = Tk. We shall further assume that the 
disk is vertically isothermal, Xfc(r, z) — Tk(r, z = 0). Models of 
dust settling sho w that most of the dust should m ostly settle within 
one scale-height dDullemond & Dom inik 2004), therefore assum- 
ing that the dust is isothermal, is at first order a reasonable assump- 
tion. The impact of this assumption will be discussed later. Under 
hydrostatic equilibrium, the resulting vertical gas distribution is a 
Gaussian 



P(r,z) = 



E(r) 



exp 



Hg{r), 



(4) 



With this definition, the gas scale height H g is: 



with k and G the Boltzmann and the gravitational constants respec- 
tively, M* the star mass, /x the mean molecular weight and mu the 
mass of the Hydrogen nuclei. H g is also a power law of the radius 

H g (r) = Ho(r/R ) h , (6) 

with the exponent h = 3/2 — q/2. The mean molecular weight \i 
is equal to 2.6 in our analysis. 

2.2 Dust Properties 

2.2. 1 Mass and grain size distributions 

Dust settling implies local changes in the dust-to-gas ratio, as well 
as local variations in the grain size distribution, whose details de- 
pend on the mechanism controlling the dust evolution. We assume 
here no radial re-distribution of dust: the dust surface density 
follows the gas surface density, and at any radius the average (i.e. 
vertically integrated) dust-to-gas ratio is equal to the standard ratio: 

X d (r)/X 9 (r) = Cstd = 1/100 (7) 

Eq|7] ensures mass conservation independently of settling. The 
value of Cstd is only a scaling factor for the total disk mass in non- 
settled disks, but also affects settling in some specific models. 

We further impose that dust settling does not change the over- 
all dust distribution as a function of grain size, and use a power law 
size distribution 



j dn(a, r, z)dz 



da 



n 



/ay 

\a J 



(8) 



(1) 



no is the number of grains at the reference size ao, dmin and 
CLmax are the minimum and maximum radius of the particles and 
Pd the exponent of the power law (usually taken from 2.5 to 4, e.g. 
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iRicci et al.ll2010bh . While the vertically integrated grain size distri- 
bution is a power law (and a fortiori, the disk averaged grain size 
distribution), because of the effect of dust settling, the local grain 
size n(a,r,z) distribution is no longer a power law of a. 

2.2.2 Dust emissivity 

The dust emissivity as a function of frequency depends on the 
dust size distribution and grain composition. Once the dust size 
distribution and grain composition are specified, several methods 
can be used to derive the emissivity values. This has to be done 
with grain sizes varying up to 5 to 6 orders of magnitude. A se- 
rious limitation is our poor knowledge of the grain composition 
and shape. Moreover, several recent observations and experiments 
show that the dust spectra l index (3 in the Far IR/mm range depends 
on the dust temperatu re dPollack et al.lll994l: lAdadze etal1ll996l : 
ICoupeaud et al.Eoil ). The Mie theory is the most popular method 
(see, e.g. lDrainell20o5) to predict dust emissivities but remains rig- 
orously exact for spherical grains only. Other methods, such as the 
Dis crete Dipole Approximation (PDA), a re heavier to handle (see, 
e.g. lDraindll988l : iDraine & Flataull20ll and references therein) 
and still suffer from the dust composition and shape limitations. 

One often uses approximate laws for the dust emissivity in 
the mm domain, such as a simple power law prescription = 
Ko (y/voY ■ Although in general applicable to the molecular clouds 
where grains remain small, in disks, this approximation is only 
valid over a relatively narrow range of frequencies. Realistic disk 
grains can result in emissivity curves which cannot be represented 
in this way at mm wavelengths, especially whe n the largest grain 
become comparable in size to the wavelength (iNatta etal.ll2004l : 
IRicci et al.|[2010bl : llsella et al.|[2009 N ) . Furthermore, in settled disks, 
such a representation would no longer be convenient, as relating 
the effective and (3 to the dust settling parameters is a non triv- 
ial task. Thus a 2-D (r,z) distribution of the dust emissivity as a 
function of wavelength needs to be computed once the settling pa- 
rameters are specified. 

Given the important unknowns in the dust geometry and com- 
position, we have elected to use a parametric method to model the 
dust emissivity as a function of grain size and wavelength. 

Our approach is based on the fact that the emissivity as a func- 
tion of frequency displays two asymptotic regimes, the small wave- 
lengths (a ^ A) where the absorption coefficient is dictated by the 
geometrical cross section, and the long wavelengths (a <C A) for 
which a power law applies. These two regimes are connected by a 
resonance region nea r A = 2na. To study the thermal structure 
of disks, llnoue et all J2009h parameterized the emissivity curves 
by only retaining the two asymptotic laws. However, at mm wave- 
lengths, the resonant region can contribute significantly to the emis- 
sivity. The detailed behaviour of this resonant region is not criti- 
cal, as integration over a size distribution will smooth out any fine 
structure: only the width and height matter. We thus elected to pa- 
rameterize the asymptotic regimes and the width and height of this 
resonant region in a simple way. The details are given in Appendix 

El 

To integrate over a given distribution in size, the distribution is 
sampled on discrete bins. We typically use two (logarithmic) bins 
per decade in size, except for the smallest sizes (below 1 /xm) where 
1 bin per decade is used because these small grains contribute very 
little to the emissivity at mm wavelengths (and are also less affected 
by settling effects). Within each bin, the size distribution is assumed 
to remain a power law with the same exponent pd as the integrated 
grain size distribution. Our selected functional for the emissivity 
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Figure 1. Contribution of various dust grains to the total emission (at 3 and 
0.5mm wavelengths), depending of their size, for an assumed size exponent 
Pd = 3. The vertical red dotted line represents approximately (depending 
of disk density) the separation between grains well mixed with the gas and 
grains starting to settle. 



k(is,cl) allows analytic integration over this truncated power law 
size distribution to derive the mean emissivity per unit mass. 

In the example presented in this article, the parameters have 
been adjusted in order to match the dust properties used by 
IRicci et all (l2010bh . The resulting emissivity per size bin are given 
in Table Q] and FigQ] shows the relative contribution of each bin to 
the total emission, for a size distribution index pd equal to 3. 



2.3 Dust Settling 

Although dust settling mechanism does not in general lead to a 
Gaussian vertical distribution of grains of a given size a, this often 
remains an acceptable approximation. Deviations from such verti- 
cal profile only occurs high above the typical scale height, i.e. in 
regions which contribute very little to the total dust mass (see e.g. 
iFromang & Nelsonll2009h . 

It is convenient to define a grain- size dependent scale height, 
Hd(a, r), and a "settling factor", s(a, r) = Hd(a, r)/H 9 {r) with 
a being the grain size. In our binned dust representation, a radius 
independent dust settling can be simulated by specifying the val- 
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Table 1. Dust Emissivity calculated from our simplified model. 



grain size 




k, (cm 2 g 


- 1 ) 




CL- 


a + 


0.5 mm 


0.8 mm 


1.3 mm 


3 mm 


0.01/iin 


30 fim 


8.33 


3.80 


1.69 


0.418 


30 fim 


100 fim 


40.1 


7.69 


1.84 


0.418 


0.1 mm 


0.3 mm 


51.4 


34.5 


11.1 


0.610 


0.3 mm 


1 mm 


8.60 


8.74 


9.52 


4.12 


1 mm 


3 mm 


2.75 


2.75 


2.58 


2.58 


3 mm 


10 mm 


0.860 


0.860 


0.826 


0.826 


10 mm 


30 mm 


0.275 


0.275 


0.270 


0.270 


30 mm 


100 mm 


0.0860 


0.0860 


0.0860 


0.0859 



ues of s n = s(a n ) for each bin. A two-bin representation (one 
layer of large grains, close to the mi d-plane, and one of sm all, 
near the disk surface) is also used by iD'Alessio et alj d2006h to 
study the impact of dust settli ng on disk SED. A small difference is 
that in ID'Alessio et ID d20Q6h the two grain categories are spatially 
separated, while in our case they would only have different scale 
heights. 

To obtain the s(a,r) value, we decided to use instead a 
more physical approximation based on the results of global 
numerical calculations d erived from theoretical approaches 
jFromang & Nelsonl 120091) which take into account ideal MRI- 
induced MHD turbulence predictions (Balbus &l?awlevi ri99L 
Il998l) as well as vertical stratification of dust and sas. 

In a Keplerian disk, the angular velocity is: 



n = Q 



-3/2 



Hi)' 

and relates to the scale height in hydrostatic equilibrium by: 



(9) 



(10) 



were C s is the (isothermal) sound speed. The dust stopping time is 
the typical time, for a particle of size a and density pd, initially at 
rest to reach the local gas velocity. In typical T Tauri protoplanetary 
disks, aerodynamic interactions between gas and solid particules 
smaller than ^10 m eters are well described by the Epstein regime 
jGaraud et al .1120041) . We have then for the dust stopping time the 
expression: 
_ pda 
Ts ~ pC s 

The main factor controlling the degree of settling is the dimension- 
less product of the dust stopping time r s by the angular velocity 
which fixes the dynamical time. When Qr <C 1, the dust particles 
are coupled to the gas. When Qr ^> 1, the dust particles are decou- 
pled from the gas and settles towards the midplane. This product is 
linked to the particle size a by: 



(11) 



Qr s (r, z) = 



27ip d a 



exp(z/H g (r)) 2 



(12) 



where the surface density T, g (r) is given by Eqs fTTl2l depending on 
which disk model is used. As this quantity is therefore inversely 
proportional to the gas surface density, in general the settling in- 
creases with radius. 

It is convenient to further approximate the effects of dust set- 
tling by relating the "settling factor" s(a, r) to the settling parame- 
ter Qt = (Qr s )(r,z = 0): 

H d (a,r) 



s(a,r) 



H 9 (r) 



/(O.ro) 



(13) 




(Qt) 



Figu re 2. Dust scale he ight as a function of (Q.r s )o (adapted from 
iFroma ng & Nelson 20(3. The black diamonds represent the values cal- 
culated by the simulations. The red and dashed blue lines indicate the two 
functions adopted in our study. 



For large grains, bubrulle et all dl995h and lCarballido etal.l d2006h 
have shown that a power law: 



H d (a,r) , _ 

s(a ' r) = -ff^r = (Jlro) 



(14) 



with a = —0-5 is a suitable function. With a similar representation, 
iPinte et ID (2008) found an exponent a = —0.05 from a multi- 
wavelength study of IM Lupi. However, their value is mostly con- 
strained by infrared data, and more specifically the silicate bands 
which are essentially sensitive to small grains. We have adopted 
the following law, which matches the previous asymptotic results 



s(a, r) 



1 



V 00c ) 



if Qto < uj c 
if Qtq > uj c 



(15) 



where uj c ^ a, the viscosity parameter, within a factor of or- 
der unity foubrulle et al.lll995h . From Fig. 2, we use in red uj c = 
6.5 10 ~ 4 , a val ue which slightly overest imates the settling effi- 
ciency found bv lFromang & Nelsonl (120091) . We will also discuss in 
Sec. l4.4.2l of the value ou c = 1.7 10 -3 , in blue, which on contrary 
tends to underestimate it. For small grains, the small difference be- 
tween ou r adopted exponen t of for small Qr and the value -0.05 
found by IPinte et alJ d2008r) is unimportant for our purpose, since 
the emission in the mm/submm domain is largely dominated by 
grains affected by dust settling, as illustrated by Figfl] 



2.4 Radiative Transfer 

We used the ray -tracer of the radiative transfer code DISKFIT 
JPietu et alj|2007h to generate brightness distributions at different 
wavelengths. As settling can only be observed at sufficiently high 
disk inclinations, special care was taken in definin g the image sam- 
pling t o limit the numerical effects, as described in lGuilloteau et alj 
(120 lib . This required to have radial and vertical cells smaller than 
0.05 AU. 
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Table 2. Disk physical parameters for a 1 Mq star. 



Physical characteristics Adopted values 



type of grains 


Moderate 3 mm) 




or Large « 10 cm) 


gas scale height 


Hydrostatic Equilibrium (Ecil5i 


Averaged Gas/Dust 


100 


Kinetic Temperature 


T k (r) = 30 (^)~ 0A Kelvin 


jjusi leiiipeidLuie 


J- dust — - 1 - k 


Reference radius 


c?^ — 1 nn att 
riQ — 1 UU AU 


Inner Radius 


Rint = 3 AU 


Inclinations 


70, 80, 85 and 90° 


Gas Surface Density (g 


.cm -2 ) Truncated disk (model 1, Eq[T} 


E*(r) 




P 


i 


Outer edge 


i? out = 100 AU 


Gas Surface Density ( 


g.cm -2 ) Viscous disk (model 2, Eq|2j 


E*(r) 


17.4 (fcPexp ((£)'-*) 


7 


0.5 


Tapered edge 


i? c = 50 AU 


Table 3. Settling factors s for the various grain size distributions 


grain size 


s(a) = H d /H 9 at Radius (AU) 


ci- a+ 


R int = 3 50 i? out = 100 


0.01 jum 10 jum 


1.00 1.000 1.000 


10 //in 30 /7,m 


1.00 0.867 0.613 


30 fim 100 /iin 


1.00 0.481 0.340 


0.1 mm 0.3 mm 


1.00 0.274 0.194 


0.3 mm 1 mm 


0.621 0.152 0.108 


1 mm 3 mm 


0.354 0.0867 0.0613 


3 mm 10 mm 


0.196 0.0481 0.0340 


10 mm 30 mm 


0.112 0.0274 0.0194 


30 mm 100 mm 


0.0621 0.0152 0.0108 



The settling factor s is calculated for pd = 3, pd = 1.5 g.cm -3 with a 
corresponding to the mean (mass weighted) grain radius, and the disk model 
described in Table [2] 

3 SIMULATIONS 

3.1 Sample of Disk Models 

The di sks parameters (T able El are representative of the disks stud- 
ied by iGuilloteau et al.1 J201lf) . The disks are in hydrostatic equi- 
librium with no vertical temperature gradient and orbit around a 
1 M star. The total (gas+dust) disk mass is 0.03 M©. 

3.1.1 Dust Settling and Emissivity 

We simulate the settling using the prescription of Eqfl3] Table [3] 
gives the corresponding settling factors and Fig|3]indicates the ap- 
parent scale height for various grain sizes as a function of radius. 

Following the formalism described in Section 2.2, dust param- 
eters we re adjusted to mimic the emissivity curves from lRicci et al.1 
see Appendix lAl The minimum grain size was 0.01 pm 
and the maximum grain size 3 mm for the moderate grain model or 
10 cm for the large grain model, with pd = 3. We took 9 grain bins 
for the moderate grains and 12 for the large ones for ensuring suf- 



i ■ ■ ' 1 ■ ■ ■ 1 ■ ■ ' r 




radius in AU 

Figure 3. Dust and gas scale heights as a function of the radius for different 
grain sizes for the settled model. The black curve also corresponds to the 
gas scale height. 



Table 4. Flux densities (mJy) of settled disks (Model 1) 



Frequency 


100 GHz 


230 GHz 


340 GHz 


670 GHz 


Moderate grains 










70° 


65 


450 


989 


3490 


80° 


48 


300 


634 


2100 


85° 


32 


180 


377 


1240 


90° 


7.6 


56 


137 


564 


Large grains 










70° 


9.6 


60 


134 


512 


80° 


9.0 


56 


123 


462 


85° 


8.1 


49 


107 


393 


90° 


2.4 


17 


39 


166 



ficient precision at the ALMA noise level. While compact minerals 
have large specific densities of pd = 3 — 4 g.cm -3 , we have chosen 
to use a smaller value pd — 1.5 g.cm -3 to account for the fact that 
(large) grains are expected to harbor a substantial ice cover and to 
be fluffy. The resulting emissivities are given in Table [T] 

3.1.2 Gas Surface Densities 

The gas surface density used to generate the settled disk model 
follows either Eqfl] (power law model, Model 1) or Eq|2] (viscous 
model, Model 2). Fig |4ll5l were obtained using the Model 1. Tables 
5 and 6 correspond to pseudo-observations using the Model 1. In 
Table 7 and 8, the pseudo observations were obtained using the 
Model 2. The resulting integrated flux densities are given in Table 

a 

3.1.3 ALMA Configuration 

The simulated brightness distributions obtained from DISKFIT 
were then processed through the regularly upgraded A LMA sim- 
ulator implemented in the GILDAS software package (iPetv et alj 
l2002h in order to produce the visibilities. 

As a first guess, we choose to simulate observations obtained 
using 50 antennas with a single antenna configuration, so that ob- 
servations at different wavelengths can be performed nearly simul- 
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Figure 4. Disks observed at 670 GHz under inclinations of 90°, 85°, 80° 
and 70°, from top to bottom. Left: settled disks. Middle: non-settled disks 
of same gas mass distribution and same amount of dust. Right: difference 
between these two models (non settled - settled). Simulations are made with 
moderate grains. The hatched ellipse is the PSF. 



taneously. A maximum baseline length of 2.3 km was used and 
the observations were assumed to be around the transit. Pseudo- 
observations of settled disks, located at declination 5 = —23°, 
have been created at four different frequencies, 100, 230, 340 and 
670 GHz (or in wavelengths: 3 mm, 1.3 mm, 0.88 mm and 0.48 mm, 
corresponding to the 4 initial ALMA bands 3,6,7 and 9). This leads 
to a spatial resolution of 0.30", 0. 13", 0.089" and 0.045" for Bands 
3,6,7 and 9, respectively. At the distance of the nearest star forming 
regions (120 - 140 pc for p Oph and Taurus- Auriga), the corre- 
sponding linear resolutions are 39-42, 16-18, 11-12 and 5-6 AU. In 
our case, we assume a distance of 140 pc. Thermal noise was added 
to the simulated uv data (corresponding to 30 min of observations 
for each frequency). The resulting image noise (point source sensi- 
tivity) are 13 /xJy at 100 GHz, 20 at 230 GHz, 30 at 340 GHz and 
111 at 670 GHz. 

Each disk has been imaged at 4 inclination angles (90°, 85°, 
80° and 70°). The resulting number of visibilities in the pseudo uv 
tables is 1096704. This number can be compared to the non reduced 
X 2 given in Tables 5 to 8. 



3.2 Prominent Effects of Settling 

To understand the effect of settling, it is useful to compare the im- 
ages of the same disk (i.e. having the same gas spatial distribution 
and mass) with or without settling. 

Figure |4]represents the expected images for disks observed at 
670 GHz, while Fig. [5] gives brightness profiles for cuts along and 
perpendicular to the disk midplane at the disk center. 

As expected, the vertical extent is smaller in the settled case, 
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Figure 5. Radial and vertical cuts in brightness temperature distribution 
(K) of disks with moderate grains, observed at 670 GHz under different 
inclinations. Left: black curves correspond to the settled model. Red curves 
correspond to the non-settled model. Right: the differences (non settled - 
settled) between these two models is shown in blue. The horizontal bar 
indicates the angular resolution. 



as well as the flaring index. At very high inclinations (only at 90° 
in our model), the r = 1 region of settled disks is reached at larger 
radial distances from the star, which are colder. This results in a 
lower brightness temperature. 

We find the same effect with large grains: their lower absorp- 
tion coefficient is partially compensated by higher column densities 
in the mid-plane due to stronger settling. The self-absorption effect 
will be smaller for less massive disks. Thus, a change in disk mass 
and a modification of the grain sizes result in different effects, in 
particular as a function of observing frequency. 

Finally, because the intrinsic aspect ratio is of order H/R < 
0.1, these opacity effects are critically dependent on the inclination 
(see Figs [4] and [5}. At 70°, the impact of dust settling becomes in 
general difficult to see. 



3.3 Inversion Process 

We investigate here potential ways to distinguish any settled disk 
from any non-settled one. Our approach is to analyze simulated 
images of settled disks with non- settled, homogeneous disk mod- 
els. Under this approach, settled disks may result in very unusual 
parameters which cannot be ascribed to "normal" non-settled disks. 
For example, the dust scale height Ho should be small, as well as 
the flaring index h, in comparison with the hydrostatic scale height. 

The resulting uv data sets were fitted by non- settled and verti- 
cally isothermal models under the assumption of power law (Model 
1, Eq[T} or exponential decay (Model 2, Eqj2j for the surface den- 
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Figure 6. Left: settled disks at 85° of inclination. Middle: results of the best 
model obtained with the non-settled disk model (Case 1). Right: difference 
between these two models (non settled - settled). Models are made with 
large grains. The hatched ellipse is the PSF. 



sity distribution. All frequencies were fitted simultaneously. Tables 
5, 6 and 8 were obtained with minimizations performed using the 
Model 1 and Table 7 using the Model 1 and Model 2). 

Non- settled disks are characterized by the following parame- 
ters: the position angle PA, the inclination i, the intrinsic param- 
eters i?int, i^out j and p (for the power law, R c , 7 for the vis- 
cous model), To,q, Ho,h and the dust characteristics. The later 
being a priori unknown, we assume the simple power law k(i/) = 
^0(^/^0)^ for the dust emissivity. We use here 1/0 = 10 12 Hz and 
«o = 0.1 cm 2 g _1 (for a dust to gas ratio of 1/100). As f3 is a free 
parameter in our analysis, the choice of 1/0 will affect k(is) at other 
frequencies, which is compensated in our analysis by adjusting the 
disk density. The derived disk density profiles E(r) (and in par- 
ticular the disk mass) is thus somewhat dependent on the assumed 
value of 1/0. 

Each pseudo-observation was fitted with 4 different non- 
settled disk models. The scale height was derived either under hy- 
drostatic equilibrium constraint or independently fitted, and dust 
emissivity exponent ft was assumed to be independent of the radius 
r, or evolving like its logarithm: 



P(r) = Pi + p r ]og(r/Ro) 



(16) 



This leads to 4 cases (see Tables [5ll6l ). Case 1 assumes hydrostatic 
equilibrium and (3 r — 0, Case 2 hydrostatic equilibrium and free 
f3 r , while Case 3 uses free scale height Ho and h with f3 r = 
and Case 4 all free parameters Ho, h and f3 r . As the impact on 
-Rint was found to be non significant in all cases, this parameter is 
ignored thereafter. The disk inclination i is recovered accurately in 
all cases (with typical error around 0.2°), but its knowledge controls 
the error bars on some critical parameters, in particular Ho and h. 
The position angle is also easily recovered, but has less influence 
than the inclination. 
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Figure 7. Radial and vertical cuts in brightness temperature distribution 
(K) at several wavelengths for disks inclined at 85° and large grains corre- 
sponding to Fig|6] Left: black curves correspond to the settled model. Red 
curves correspond to the non-settled model. Right: the differences (non set- 
tled - settled) between these two models is shown in blue. The horizontal 
bar indicates the angular resolution. 



4 DISCUSSION 

4.1 Analysis of the Inversion Process 

Tables l5ll8l show the results of the inversion process. In Tables|5]and 
[6] both pseudo-observations and models for the minimizations use 
the truncated disk surface density (Model 1), with grains of moder- 
ate size in Tableland large grains in Table |6] In Table IT] pseudo- 
observations, made with the viscous Model 2 and containing large 
grains, are analysed by both Models 1 and 2. Finally, Table prefers 
to pseudo-observations obtained with Model 1 and fitted using the 
Model 2, for moderate size grains. 

The case with grains of moderate size illustrates best the prob- 
lems. It leads to rather strong continuum flux (Table |4j, and the op- 
tically thick zone is sufficiently large to measure directly the dust 
temperature from the surface brightness. The formal errors are very 
small, indicating that thermal noise is not a limitation here. 



4.1.1 Deriving the Scale Height 

When viewed edge-on, the hydrostatic equilibrium assumption 
(cases 1 and 2) leads to unusual results. The derived temperature 
is forced towards low values to better mimic the small disk thick- 
ness (rsj 19 K instead of 30 K). This is also true when minimizing a 
Model 1 by a Model 2 (Table [8]). A side effect is an apparent radial 
dependency of the dust emissivity index (f3 r / 0) which is due to 
the nonlinearity of the Planck function at low temperatures. Relax- 
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Table 5. Minimizations of a settled disk (model 1) by an homogeneous one (model 1): moderate size grains 



Disk Case & 


T 


<2 


V 


^out 


Pi 


Pr 


^0 




2 

X 








1.0 


100 AU 






14.6 AU 


1.30 




(Case 1) 




















70° 


30.3 ± 0.08 


0.396 ± 0.001 


0.966 ± 0.002 


95.8 ± 0.02 


0.613 ± 0.001 


10] 


(14.7) 


(1.30) 


1458941 




33.8 ± 0.07 


0.378 ± 0.002 


1.32 ± 0.004 


96.1 ± 0.03 


0.420 ± 0.001 


10] 


(15.5) 


(1.31) 


1696971 




26.6 ± 0.10 


0.629 ± 0.007 


0.692 ± 0.006 


97.9 ± 0.04 


0.341 ± 0.002 


10] 


(13.7) 


(1.18) 


1661490 


90° 


18.8 ± 0.22 


0.810 ± 0.035 


-1.20 ± 0.03 


101.0 ± 0.07 


0.892 ± 0.005 


10] 


(11.6) 


(1.10) 


1332034 


(Case 2) 




















70° 


30.3 ± 0.08 


0.396 ± 0.001 


0.965 ± 0.006 


95.8 ± 0.02 


0.612 ± 0.002 


0.000 ± 0.005 


(14.7) 


(1.30) 


1458923 




33.8 ± 0.13 


0.378 ± 0.001 


1.32 ± 0.01 


96.1 ± 0.03 


0.420 ± 0.003 


0.003 ± 0.008 


(15.5) 


(1.31) 


1697294 




26.4 ± 0.17 


0.640 ± 0.005 


1.16 ± 0.01 


98.4 ± 0.06 


0.230 ± 0.003 


-0.23 ± 0.01 


(13.7) 


(1.18) 


1660374 


90° 


19.3 ± 0.3 


0.622 ± 0.013 


0.55 ± 0.07 


102 ± 0.10 


0.440 ± 0.005 


-1.21 ± 0.02 


(11.7) 


(1.19) 


1325182 


(Case 3) 




















70° 


28.3 ± 0.04 


0.423 ± 0.001 


1.09 ± 0.005 


99.7 ± 0.007 


0.686 ± 0.001 


10] 


2.55 ± 0.05 


1.06 ± 0.016 


1205325 


80° 


29.2 ± 0.03 


0.412 ± 0.001 


1.43 ± 0.009 


99.7 ± 0.009 


0.704 ± 0.002 


10] 


2.49 ± 0.04 


0.94 ± 0.007 


1193405 




30.2 ± 0.03 


0.391 ± 0.001 


0.27 ± 0.06 


100. ± 0.03 


0.682 ± 0.007 


10] 


2.91 ± 0.06 


1.18 ± 0.03 


1155545 


90° 


31.2 ± 0.19 


0.39 ± 0.034 


-1.20 ± 0.05 


101. ± 0.07 


0.833 ± 0.007 


10] 


3.06 ± 0.05 


-0.09 ± 0.02 


1117094 



28.1 ± 0.04 

29.1 ± 0.03 

30.2 ± 0.03 
30.2 ±0.19 



0.425 ± 0.001 
0.415 ± 0.001 
0.391 ± 0.001 
0.54 ± 0.038 



1.22 ± 0.008 
1.78 ± 0.02 
0.35 ± 0.09 
-1.65 ± 0.09 



99.8 ± 0.007 
99.8 ± 0.008 

100 ± 0.06 

101 ± 0.07 



0.648 ± 0.002 
0.644 ± 0.003 
0.601 ± 0.007 
0.897 ± 0.006 



-0.080 ± 0.005 
-0.21 ± 0.02 
-0.46 ± 0.07 
0.17 ± 0.02 



2.49 ± 0.06 
2.52 ± 0.04 
2.95 ± 0.07 
3.05 ± 0.06 



1.05 ± 0.014 
0.97 ± 0.007 
1.26 ± 0.03 
-0.11 ± 0.03 



1205325 
1192868 
1155401 
1117059 



Numbers between brackets [] indicate fixed parameters. Numbers between parentheses are derived from another parameter ( Hp from T^, and h from q under the hydi 

the parameters. See section |33l for the definition of Cases. 



equilibrium hypothesis). The second row indicates the expected values of 



Table 6. Minimizations of a settled disk (model 1) by an homogeneous one (model 1): large grains 



Disk Case & 


T 


Q 


P 


^out 


Pi 


Pr 


^0 


h 


x 2 


inclination 


30 K 


0.4 


1.0 


100 AU 







14.6 AU 


1.30 




(Case 1) 




















70° 


22.2 ± 1.2 


0.453 ± 0.014 


0.92 ± 0.013 


99.1 ± 0.2 


0.337 ± 0.004 


10] 


(12.6) 


(1.27) 


1099461 


80° 


19.2 ± 0.6 


0.555 ± 0.011 


0.78 ± 0.01 


96.5 ± 0.2 


0.366 ± 0.004 


10] 


(11.7) 


(1.23) 


1103372 




17.8 ± 0.3 


0.668 ± 0.010 


0.70 ± 0.01 


95.8 ± 0.2 


0.464 ± 0.005 


10] 


(11.2) 


(1.17) 


1115798 


90° 


11.9 ± 0.4 


0.967 ± 0.055 


-0.82 ± 0.06 


96.6 ± 0.3 


0.98 ± 0.02 


10] 


(9.1) 


(1.02) 


1113205 


(Case 2) 




















70° 


15.0 ± 1.0 


1.10 ± 0.04 


-0.07 ± 0.04 


99.8 ± 0.1 


0.45 ± 0.008 


0.126 ± 0.004 


(10.3) 


(0.95) 


1099558 




13.8 ± 0.9 


1.30 ± 0.017 


-0.45 ± 0.02 


97.9 ± 0.2 


0.50 ± 0.007 


0.182 ± 0.004 


(9.9) 


(0.85) 


1102652 




17.7 ± 0.24 


0.674 ± 0.010 


0.73 ± 0.02 


95.8 ± 0.2 


0.45 ± 0.008 


-0.01 ± 0.01 


(11.2) 


(1.17) 


1114358 


90° 


13.0 ± 0.35 


0.660 ± 0.027 


0.73 ± 0.08 


97.8 ± 0.2 


0.43 ± 0.02 


-0.93 ± 0.05 


(9.6) 


(1.17) 


1112144 


(Case 3) 




















70° 


32.0 ± 1.0 


0.40 ± 0.02 


1.01 ± 0.02 


99.6 ± 0.06 


0.266 ± 0.004 


10] 


1.6 ± 0.7 


1.43 ± 0.13 


1097745 




33.0 ± 0.9 


0.35 ± 0.02 


1.08 ± 0.02 


98.9 ± 0.2 


0.270 ± 0.004 


10] 


1.9 ± 0.5 


0.68 ± 0.07 


1098468 


85° 


33.9 ± 1.0 


0.32 ± 0.02 


1.13 ± 0.02 


98.7 ± 0.2 


0.268 ± 0.004 


10] 


1.4 ± 0.2 


0.53 ± 0.12 


1099739 


90° 


22.3 ± 0.8 


1.49 ± 0.07 


-1.34 ± 0.07 


100. ± 0.2 


0.51 ± 0.02 


10] 


1.4 ± 0.2 


0.03 ± 0.07 


1096614 


(Case 4) 




















70° 


27.8 ± 1.3 


0.510 ± 0.05 


0.84 ± 0.08 


99.7 ± 0.06 


0.30 ± 0.01 


0.027 ± 0.006 


1.1 ± 0.3 


1.3 ± 0.4 


1097740 




35.9 ± 1.3 


0.310 ± 0.03 


1.17 ± 0.03 


98.9 ± 0.14 


0.23 ± 0.01 


-0.028 ± 0.005 


1.6 ± 0.6 


0.6 ± 0.2 


1098426 




34.5 ± 1.1 


0.309 ± 0.03 


1.20 ± 0.03 


98.8 ± 0.2 


0.22 ± 0.01 


-0.043 ± 0.007 


1.3 ± 0.2 


0.45 ± 0.2 


1099695 


90° 


58.0 ± 3.1 


-0.28 ± 0.03 


1.64 ± 0.08 


100. ± 0.2 


0.11 ± 0.01 


-0.64 ± 0.05 


1.3 ± 0.2 


0.05 ± 0.1 


1096561 



ing the hydrostatic equilibrium hypothesis (cases 3 and 4) allows 
us to recover the input temperature profile. 

Tables \5\ [6] and [8] also show that the constraint from the ap- 
parent thickness is less important than that from the dust temper- 
ature, so that the fitted scale height in the hydrostatic equilibrium 
hypothesis remains unduly large (Fig.[6]and FigJTJ. This artificially 
creates a deficit of emission close to the mid-plane and an excess 
at high altitude. For cases 3 and 4, there is a lack of flaring at all 
at 90°: the settled disk is best fitted with a constant thickness. At 
less extreme inclinations, however, disks appear mildly flared. Not 
only Ho is constrained, but the apparent flaring index h also devi- 
ates quite significantly in the settled case from the initial value (1.3 
in our model, a range between 1.1 - 1.5 being expected for most 
disks). 



4.1.2 Spectral Index 

Our settled disk are composed of several populations of grains. 
Each grain population has its own spectral index f3. If the whole 
dust emission was optically thin and homogeneously distributed, 
a mean f3 (defined as the spectral index between two wavelengths 
only: 0.5 and 3 mm) of ~ 0.61 for the moderate grains (Table [5} 



and ~ 0.29 for the large grains (Table [6} is expected from the opac- 
ity curves in Appendix lAl The fitted f3 is often different because /3 
is not an intrinsic parameter of the dust: our assumed dust proper- 
ties cannot be represented by a single power law between 3 and 0.5 
mm, but exhibit a more complex behaviour (see Appendix). The fit- 
ted (3 is more affected for edge-on disks, because the flux densities 
at each frequency strongly depend on the degree of settling, thus 
affecting the relative weights of each observation. 



4.1.3 Degeneracy between f3 r and p 

At very high inclinations (e.g. 90°), settling increases the opacity 
in the disk plane. A fit of a constant f3 (f3 r = 0) leads to a value 
of the exponent p of the radial density profile driven towards neg- 
ative values, to offer sufficient self-absorption from the cold outer 
regions. The independent fit of Ho and its exponent h (case 4) is 
not sufficient to compensate this effect. Although this suggests that 
viscous-like surface density profiles (see Eq|2j with negative 7 may 
better fit the images, this is not the case because such profiles drop 
too sharply after their critical radius R c . Furthermore, there is some 
"hidden" degeneracy between p and f3 r and the minimization pro- 
cess may converge towards one or the other solution. 
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Table 7. Tapered Edge disk: large grains 



Disk Case & 
inclination 


^0 
30 K 


<? 

0.4 


P 


-Rout 


Pi 


P r 



^0 
14.6 AU 


h 

1.30 


x 2 


70° 
(Case 1) 
(Case 2) 
(Case 3) 
(Case 4) 


32.9 ± 0.1 
29.8 ± 0.07 
28.7 ± 0.11 
25.4 ± 0.07 


0.310 ± 0.001 
0.348 ± 0.001 
0.379 ± 0.002 
0.426 ± 0.001 


1.73 ± 0.002 
2.16 ± 0.004 
1.65 ± 0.002 
1.90 ± 0.006 


91.7 ± 0.07 
92.9 ± 0.07 
94.1 ± 0.08 
94.0 ± 0.08 


0.689 ± 0.002 
0.419 ± 0.004 
0.725 ± 0.002 
0.629 ± 0.005 


[0] 

-0.322 ± 0.004 
[0] 

-0.181 ± 0.004 


(15.3) 

(14.6) 
20.0 ± 0.13 
18.4 ± 0.14 


(1.35) 

(1.33) 
2.32 ± 0.007 
2.31 ± 0.009 


1151874 
1141912 
1143237 
1134115 






As Table^r ; 


i viscous (model 2) settled disk fitted by an homogenous disk with sharp edge (model 1). 






Disk Case & 
inclination 


30 AU 


Q 

0.4 


P 
0.5 


Rc 
50 AU 


Pi 


Pr 



^0 
14.6 AU 


h 

1.30 


x 2 


70° 

(Case 1) 
(Case 2) 
(Case 3) 
(Case 4) 


30.5 ± 0.09 

29.6 ± 0.08 
25.1 ± 0.10 
25.3 ± 0.09 


0.387 ± 0.001 
0.391 ± 0.001 
0.473 ± 0.002 
0.463 ± 0.002 


0.479 ± 0.003 
0.578 ± 0.003 
0.459 ± 0.003 
0.520 ± 0.003 


49.3 ± 0.04 
46.2 ± 0.09 
51.2 ± 0.05 
49.2 ± 0.09 


0.620 ± 0.002 
0.536 ± 0.003 
0.682 ± 0.002 
0.628 ± 0.004 


[0] 

-0.092 ± 0.002 
[0] 

-0.052 ± 0.002 


(14.7) 
(14.5) 
3.4 ± 0.3 
2.9 ± 0.4 


(1.31) 

(1.30) 
1.73 ± 0.03 
1.65 ± 0.04 


1106571 
1106488 
1099314 
1099013 



As Table[5Tfor a viscous (model 2) settled disk fitted by an homogenous viscous disk (model 2). 



Table 8. Tapered Edge disk: moderate size grains 



Disk Case & 


T 


<? 


P 


^out 


Pi 


Pr 


^0 


h 


x 2 


inclination 


30 


0.4 






0.613 





14.6 


1.30 




(Case 1) 




















70° 


34.8 ± 0.07 


0.319 ± 0.001 


2.14 ± 0.002 


103.0 ± 0.03 


0.442 ± 0.001 


10] 


(15.7) 


(1.34) 


1276606 


80° 


35.0 ± 0.07 


0.346 ± 0.002 


1.85 ± 0.004 


103.0 ± 0.03 


0.292 ± 0.001 


10] 


(15.8) 


(1.33) 


1415333 




29.5 ± 0.13 


0.504 ± 0.004 


1.31 ± 0.006 


109.0 ± 0.06 


0.293 ± 0.003 


10] 


(14.5) 


(1.25) 


1683039 


90° 


19.9 ± 0.26 


0.542 ± 0.008 


0.16 ± 0.01 


109.6 ± 0.09 


0.930 ± 0.006 


10] 


(11.9) 


(1.23) 


1319208 


(Case 2) 




















70° 


35.5 ± 0.07 


0.306 ± 0.003 


2.05 ± 0.005 


108 ± 0.03 


0.484 ± 0.005 


0.156 ± 0.003 


(15.9) 


(1.35) 


1299320 ** 




34.8 ± 0.07 


0.349 ± 0.004 


1.70 ± 0.007 


111 ± 0.03 


0.483 ± 0.008 


0.300 ± 0.005 


(15.7) 


(1.33) 


1393296 


85° 


28.6 ± 0.09 


0.533 ± 0.005 


1.35 ± 0.008 


116 ± 0.05 


0.334 ± 0.010 


-0.016 ± 0.007 


(14.3) 


(1.23) 


1679815 


90° 


20.6 ± 0.15 


0.550 ± 0.009 


0.18 ± 0.02 


1 1 1 ± 0.07 


0.398 ± 0.021 


-1.25 ± 0.013 


(12.1) 


(1.22) 


1312764 


(Case 3) 




















70° 


35.8 ± 0.1 


0.306 ± 0.001 


2.23 ± 0.002 


106 ± 0.03 


0.466 ± 0.001 


10] 


13.4 ± 0.04 


1.33 ± 0.005 


1262191 




35.0 ± 0.06 


0.302 ± 0.002 


2.26 ± 0.005 


110 ± 0.03 


0.365 ± 0.001 


10] 


6.6 ± 0.05 


1.19 ± 0.006 


1187642 


85° 


32.8 ± 0.08 


0.320 ± 0.003 


2.69 ± 0.009 


122 ± 0.05 


0.461 ± 0.002 


10] 


2.8 ± 0.04 


0.53 ± 0.006 


1129793 


90° 


29.2 ± 0.13 


0.294 ± 0.012 


0.20 ± 0.01 


113 ± 0.06 


0.754 ± 0.076 


10] 


4.0 ± 0.04 


0.10 ± 0.02 


1162383 


(Case 4) 




















70° 


35.9 ± 0.06 


0.305 ± 0.002 


2.14 ± 0.004 


108 ± 0.04 


0.493 ± 0.005 


0.099 ± 0.003 


13.3 ± 0.03 


1.32 ± 0.008 


1278814** 




33.5 ± 0.08 


0.330 ± 0.003 


2.25 ± 0.010 


114 ± 0.04 


0.372 ± 0.005 


0.172 ± 0.007 


5.57 ± 0.05 


1.10 ± 0.005 


1186973 


85° 


32.4 ± 0.09 


0.327 ± 0.005 


2.68 ± 0.017 


127 ± 0.06 


0.451 ± 0.010 


0.268 ± 0.02 


2.35 ± 0.03 


0.44 ± 0.007 


1124899 


90° 


34.4 ± 0.10 


0.068 ± 0.004 


1.21 ± 0.007 


122 ± 0.05 


0.738 ± 0.008 


-0.380 ± 0.005 


4.04 ± 0.03 


0.34 ± 0.016 


1145282 



As Tablef5lfor a viscous (model 2) settled disk fitted by an homogenous disk with sharp edge (model 1). ** Results probably not converged, as their %2 is greater than that of the simpler (3 r = case. 



4.1.4 Impact of the Surface Density Profile 



4.1.5 Consequences 



Table |5] suggests that the scale height can be apparently con- 
strained independently of the temperature profile even at moder- 
ate (i = 70°) inclination (basically all input parameters are recov- 
ered properly). This result is due to the assumed sharp truncation at 
i^out = 100 AU (Model 1). The apparent (projected) width of this 
sharp edge is a strong indicator of the actual disk thickness. 

Table |7] shows results for pseudo-observations obtained with 
a more realistic continuous profile (Model 2, with 7 = 0.5 and 
R c = 50 AU). When fitted by a Model 1 (top panel of Table [7]), 
the required scale height is large and the flaring index reaches non 
physical values of the order of 2.5. This is an attempt to fit the emis- 
sion beyond the derived outer radius. On the contrary when fitted by 
a Model 2 (bottom panel), a small scale height is indeed recovered. 
This result indicates that at inclinations below 80°, the recovered 
scale height is sensitive to the exact shape of the surface density 
distribution, and cannot in general be determined accurately. Ta- 
ble [8] shows results of tapered disks (Model 2) fitted by a truncated 
power law for different inclinations (Model 1) and moderate size 
grains. At inclinations > 80°, the differences between the true disk 
density structure and the one assumed in the analysis do not signif- 
icantly affect the derivation of the scale height. Other parameters, 
such as the temperature, are somewhat affected by the improper 
surface density profile rather than by settling. 



For the large grains models, the flux densities and the optical depths 
are lower. The same trends are found. Large grains settle more effi- 
ciently and the fitted scale height is even smaller than in the previ- 
ous case. Since the optically thick core is small, some degeneracies 
start appearing between T, q and as a purely optically thin 
emission only depends on ST and p + q. 

In all cases, the inconsistencies appearing when fitting by a 
standard, non-settled disk model, clearly flag the "observed" disk 
as being unusual, and combined with the low absolute values of the 
scale height (~ 2 — 3 AU), point towards dust settling as the only 
reasonable cause of the discrepancies. Moreover, settled disks ac- 
tually appear "pinched" (h < 1) rather than flared (h > 1). The 
above analysis also demonstrates that radius dependent settling as 
derived from MRI simulations and theoretical analysis can be dis- 
tinguished from radius independent one, the later would not affect 
the flaring index value h. However, directly retrieving the settling 
factor s(a,r) will remain largely model dependent, as this would 
imply to deconvolve from the grain size distribution n(a), which 
remains unknown. Even with prior knowledge of n(a), the strong 
smoothing resulting from this size distribution would severely limit 
the capability to retrieve s(a, r) from H(r). 

Others parameters like ft, (3 r or p are sensitive to the dust set- 
tling at inclinations ^ 80° but can only serve as secondary indica- 
tors. In real data, the f3 r which deviates from its original value at 
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inclinations > 80° may be due either to dust settling o r to radial 
variations of the grain properties (Guilloteau et al 

An inclination close to 90° is clearly the more suitable case to 
study settling since opacity and brightness temperature effects are 
maximum. Taken into account the various uncertainties, in partic- 
ular on the surface density and radial grain properties, our results 
suggest that observations of settling would be possible at inclina- 
tions > 75 - 80°. 



4.2 Impact of the various Wavelengths 

The above studies show that all the impact of dust settling is only 
in the effective scale height (Fig (8} and a priori we may expect 
that the highest frequency data, which has the highest spatial res- 
olution, may be sufficient in itself. This must be moderated by 
a number of caveats, however. First, the best signal to noise de- 
pends on the dust properties and is not necessarily at the highest 
frequency. Second, the apparent (geometrically constrained) scale 
height must be compared to the hydrostatic scale height to prove 
settling. This implies that a) the (gas) temperature or the dust tem- 
perature as a proxy, should be known, and b) the stellar mass must 
also be constrained to a reasonable accuracy (to derive H^, see 
Eq|5]). In principle, the gas temperature can be retrieved by imag- 
ing thermalized lines. However, as most chemical models predict 
that simple molecules lies in a layer about 1-2 scale height above 
the disk plane (becau se of depletion on dust grains in the cold 
denser regions, see e.g. lSemenov & Wiebel201ll) . finding a suitable 
probe for the disk plane is not straightforward. In our approach, the 
dust temperature is derived by resolving the optically thick parts 
of the dis k. With radial gradients of the dust emissivity index like 
found by iGuilloteau et al J fcOllh and predicted by simulations of 
iBirnstiel et alJbOlOh . the proper identification of an optically thick 
core region requires at least 3 frequencies. Thus unless some gas 
temperature can be derived independently, a 3-wavelength study 
seems required to avoid ambiguities in identifying dust settling. 

The relative ability of each of our 4 observing wavelengths can 
be evaluated. For the two shortest ones (0.5 and 0.8 mm), the errors 
on the derived parameters (e.g. Ho and h) approximately scale as 
the wavelength. Since the signal-to-noise ratio is similar at both 
frequencies, the driving factor is the angular resolution. The errors 
then strongly increases for 1.3 mm, which no longer has sufficient 
resolution, while the 3 mm data are practically unable to provide 
any quantitative constraint. Good observing conditions at 0.8 mm 
data being much more frequent than at 0.5 mm, this wavelength 
may be the best compromise in term of sensitivity and angular res- 
olution if only one wavelength can be observed. 

We note that the error on To in the combined analysis is lower 
than the simple weighted average of the 4 independent determina- 
tions which shows the gain in the multi- wavelength approach. 

We finally made a last check at 3 mm on long baselines by 
using the Model 1 to produce pseudo-observations in the case of the 
moderate grain size distribution and assuming an inclination angle 
of 85°. Baseline lengths of ~ 11 km provide an angular resolution 
of about 0.06" or 8 AU, similar to that reached at 0.5 mm with the ~ 
2 km baselines. We mimic 4 hours of observations. We analyzed the 
pseudo-observations using the non settled disk Model 1 and found 
that the scale height is marginally constrained with Ho = 2.2 ± 0.7 
AU and h = 1.16 ± 0.3. This also barely differs (by ~ 1.5a) 
from the 4 wavelength fit where we obtain Ho = 3.13 ± 0.03 AU. 
The large errors at 3 mm are due to insufficient sensitivity. Thus, 
measuring the differential settling between 0.5 and 3 mm would be 
very time consuming. 




70 80 90 

Disk inclination 

Figure 8. Derived apparent scale heights and flaring index as a function of 
disk inclination. All results show viscous pseudo-observations (Model 2) 
fitted by truncated disks (Model 1), using grains of moderate size, with free 
scale height and radial dependent ft (case 4). Stars show results obtained 
with the strong dust settling prescription while filled symbol are obtained 
with the lowest one. The dashed lines indicate the expected hydrostatic scale 
height and flaring index for non-settled disks. 



4.3 Comparison with other imaging simulations 

Using the code MC3D, ISauter & Wol3 (1201 ll) have investigated 
dust settling by producing intensity maps of dust disks from 1.0 /im 
up to 1.3 mm. Their analysis differs from ours in three major points. 

First, they only assume two dust grain distrib utions (small and 
large) following the parametrization proposed by iD'Alessio et al.l 
(2006). Their small grain population is ISM-like, the large grain 
distribution extends up to a ma x = 1 mm. This parametrization is 
similar to the 2-bin version of our radius independent settling 
models (Section EJ). It is very well suited to study settling in the 
NIR and Mid-IR because of the high dust opacity but has a too 
small number of bins to properly mimic dust settling at mm wave- 
lengths. The maximum grain size may not be sufficient, as shown 
for example in Fig(3]where the larger grains significantly contribute 
to the mm emissions. They also only use stronger settling parame- 
ters, with their large grain scale height smaller than the small grain 
one by factors 8, 10 or 12. This roughly corresponds to the settling 
factor in our large grain case, but the ratio is of the order of 3 for 
our less extreme grain sizes. 

Second, they do not take into account the ALMA transfer 
function. This is adequate only with sufficient uv coverage, which 
is not obtained with short integrations on very long baselines. 

Third, and most importantly, they only compare the settled 
model with the non- settled disk in four positions. Such a method, 
optimized for IR data, does not use all the information contained in 
the maps or ALMA observations. Moreover, as the dust opacity is 
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changing with the wavelength, the optimum positions should vary 
accordingly. 

Given these differences, comparisons are not straightforward. 
As expected, we both find weaker flux and reduced flaring for 
highly inclined settled disks, but our method appears much more 
discriminant and applicable to a wider range of disk inclinations. 

4.4 Critical discussion 

4.4.1 Temperature Structure: 

We assume that th e temp erature is vertically isothermal (as did 
iFromang & Nelson] J2009h ). In real disks, the temperature is ex- 
pected to rise two or three scale heights above the disk plane. Dust 
settling will affect this temperature gradient which is mostly driven 
by the distribution of small to mid-size grains (~ 0.1 to 10 /xm) be- 
cause they control the opacity to incoming radiation. These grains 
exhibit only limited settling. Indeed, because the apparent scale 
height at mm wavelengths is a factor 3 to 4 times smaller than the 
hydrostatic scale height, more than 99 % of the mm flux is built 
in within one hydrostatic scale height, in which temperature gra- 
dients should be negligible. The location of the super-heated layer 
changes with dust settling, but not to the point where it will sub- 
stantially 50 %) affect the temperature within one pressure scale 
height 

lHasegawa & Pudritzl j201lh have recently stu died the effe ct of 
dust settling on the dust temperature using MC3D (Wolf 2003). At 
5 AU from the star, they found that the dust temperature near the 
mid-plane (within z ^ 0.3 scale-height) is somewhat lower in a 
settled disk than in a well-mixed one (see their Fig. 4). The super- 
heated layer appears however hotter (60 K instead of 40 K). The 
thickness of the impacted cold layer is around 0.5 AU, much too 
small even for the longest ALMA baselines. At the larger radii (50 
to 100 AU) investigated in our study, the impact of dust settling on 
the temperature structure will be much less significant, because the 
temperature gradients scale with the dust opacity (oc 1/r in typical 
disks), as well as with incoming radiation flux (oc 1/r 2 ). 

Vertical temperature gradients are however expected to play 
a role in the apparent s cale height at optica l or NIR wave- 
lengths. Indeed for HH30, iBurrows et al.l dl996h derived a much 
larger scale height from 2 /xm scattered light using the HST than 
IGuilloteau etafl J2008h from IRAM PdBI data: this is more likely 
a manifestation of temperature gradient than of dust settling. 



4.4.2 Settling Shape and Viscosity Parameter: 

We have tested a prescription of the settling wh ich has been derived 
from M RI driven turbulence simulations from IFromang & Nelson] 
(120091) . These simulations span a limited range of (Qr s )o, and Fig- 
ure |2 suggests that other settling factors may be used. We also per- 
formed simulations with u c about twice larger (dashed blue curve 
on Fig. 2), leading to smaller settling but without any major change 
in our results as can be seen on Fig. 8. The dust scale heights are 
affected by at most 30-40 %, but are still strongly smaller than 
the gas (~ 3 AU instead of 15 at R = 100 AU for the grain range 
size 0. 1-lmm) and still easily distinguished from the unsettled case. 
Furthermore, the settling degree is also directed linked to the dust 
specific density, which is generally assumed to be between 1 and 
3 g.cm -3 . As our grains have a relatively low dust specific den- 
sity (1.5 g.cm -3 ), and then are more coupled to the gas, the dust 
settling degree we generally used can be considered as medium. 
The measurable effects on the apparent flaring index indicate that 



the settling produced by MRI can be distinguished, in some cases, 
from a radially constant settlin g. For instance, in the simulation 
from dFromang & Ne lson 2009), there is no dead zone, leading to 
an underestimate of the dust settling in the inner disk (r < 10 AU). 

Having measured Hd/H g , it is tempting to directly quan- 
tify the viscosity parameter a. When this ratio is inferior to 1 , 
the settling efficiency is related to it by dDubrulle et all Il995l : 
ICarballido et 



iinciency 

"all kood) : 



s(a,r) 



H g {r) 



(17) 



However, we do not measure Hd(a,r)/ H g (r) as a function of 
grain size, but only an ensemble averaged with an a priori unknown 
size distribution, and a weighting function depending on the dust 
emissivity as function of size and wavelength. Furthermore, Qr 
scales as the inverse of the gas surface density, which varies by fac- 
tor of a few across the disk radius. These two effects are difficult 
to separate from the direct impact of a in the above formula. Thus, 
Hd/Hg strongly depends on many parameters and quantifying a 
in this way appears impracticable. This conclusion is unfortunately 
re-inforced by the large integration times which would be required 
to measure differential settling between 0.5 and 3 mm, a minimal 
step to attempt any correction from the grain size distribution. 

Each grain size bin is assumed to follow a Gaussian shape 
vertical d istribution. Deviation from Gaussianity are expected, as 
shown bv IFromang & Nelson! (120091) from their MRI simulations. 
However, as these deviations occur above two-three scale-heights, 
they play a minor role in the mm/submm results, like the tempera- 
ture profile. 

4.4.3 Disk Size: 

We use a disk outer radius of 100 AU, or a similar characteristic size 
for the tapered-edge profile. This is consistent wit h the sizes found 
for disks in the Taurus Auriga or p Oph regions (IGuilloteau et al.l 



2QlTI: [Andrews et al The median disk outer radius in 



Guilloteau etail (1201 lh is 13 AU 
alj (|2 



Guilloteau et al. M2011I) also show that large grains are essen- 



tially located in the inner 70 to 100 AU. The existence of a ra- 
dial gradient in the grain size distribution will affect the radial de- 
pendency of settling. The presence of smaller gr ains beyond 100 
AU, a s suggested by the observational results of (iGuilloteau et al.l 
l201hV will increase the apparent scale height there. These outer re- 
gions contribute to less than 10 to 30 % of the total flux, so this 
increased scale height will not mask the settling from the inner re- 
gions. If grain growth is maximal near the star, the settling will 
become more important in the inner regions. Thus, the radial gra- 
dient of grain size should increase the expected values of the flar- 
ing exponent h. While this reduces one of the signature of settling, 
the effect can easily be mitigated by a slightly modified analysis 
method, as the change in h is correlated with the change in grain 
size. This is amenable to simple parametrization, at the expense of 
one additional parameter. 

4.4.4 Disk Mass: 

We adopted a disk mass of Md = 0.03 M which corresponds to 
flux densities given in Table |4] and are similar to those of disks 
found in e.g. the Taurus region (a factor 2 larger for the small 
grain case and a factor 2 smaller for the large grain case). With 
enough sensitivity, it is preferable to observe disks with low fluxes 
for two reasons: 1) at mm wavelength, low flux densities can be 
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a sign of larger grains (a > A) which settle more efficiently and 
2) low flux densities may indicate lower densities which decrease 
the coupling between gas and dust. In the "large grain" approxi- 
mation, the settling factor s scales as \J Md/a (from Eqflj] and 
Q3}, while for optically thin emission, the flux density S u scales as 
K,(a)Md ~ Md/a, so that s oc \/S^, the settling efficiency 1/s 
may actually be larger for disks with lower observed flux densities, 
unless these low flux values are just due to disks of similar intrinsic 
densities, but smaller outer radii. 

4. 4. 5 Disk Radial Structure: 

We performed a series of tests where we fit the settled model (using 
moderate and large grain size distributions) with a viscous law for 
the surface density (Model 2) by a non- settled model assuming a 
power law surface density (Model 1). The results of the minimiza- 
tion (Table[8j show that even if the limited knowledge of the density 
profile unavoidably affects the precision with which settling is con- 
strained, it does not mask its existence. For disks inclined by more 
than 80° , the derived scale height and its flaring clearly exhibit the 
behavior expected in case of dust settling. This result is strongly 
encouraging, especially as the radial density profile of disks is still 
a debated issue. 



4. 4. 6 Disk Geometry: 

Small departures from perfect geometry and rotational symmetry, 
like warps and spiral patterns, may affect our ability to constrain 
the scale he ight. Mis-alignm ents between jets and disks by a 1-2° 
(e.g., HH 3 Q|Petv et alJl2006h. a nd warps of similar magnitude (e.g. 
P Pictoris jMouillet et all 19971) are known to exists. They may ulti- 
mately limit the apparent scale height to about H/R~ 0.03, which 
is comparable to our "normal" grain size case. The very strong set- 
tling predicted for large grains by the MRI turbulent model may be 
beyond reach because of this practical limitation. 



4.4.7 Instrumental Effects: 

We have shown that thermal noise is not a major limitation to mea- 
sure dust settling. We investigated here the impact of atmospheric 
phase noise by adding antenna based, Gaussian distributed phase 
errors, as would be expected after radiometric phase correction. 
Figure [9] shows the impact on the results. As expected, the impact 
is much worse on the brightest sources (the small grain case). How- 
ever, for reasonable observing conditions (antenna based rms noise 
below 30°), phase noise does not prevent the measurement of the 
scale height. The flaring index is more affected, but h still shows 
significant deviations at 85° of inclination. 



5 SUMMARY 

We have studied how ALMA can be used to quantify the degree 
of dust settling in proto-planetary disks around T Tauri stars. We 
simulated settled disks using prescriptions for dust settling based on 
MRI driven viscosity. Using a parametric model to fit the predicted 
dust emission as a function of wavelength, we show to what extent 
settling can be constrained. Our main findings are 

• For the characteristic dust disk sizes found by previous mm 
surveys, dust settling can be measured in typical disks with moder- 
ate integration times (of about 2 hours per source), using baselines 
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Figure 9. Impact of the phase noise on the measurements of settling. Top: 
reduced x 2 , middle: flaring index h, bottom: apparent scale height at 100 
AU. Left panels are for the small grain case, and right panels for the big 
grain case. The truncated power law (model 1) is used. Crosses are for i = 
70°, and squares for i = 85°. The horizontal lines indicate the expected 
values for un-settled disks, derived from the hydrostatic conditions. 



of the order of 2 km at the distance of the nearest star forming re- 
gions (120 - 140 pc). 

• This is possible only for disks more inclined than ~ 75 — 80° . 

• Unless the gas scale height can be independently derived, at 
least 3 frequencies are needed to unambiguously identify settling, 
by comparing the apparent scale height to the derived dust temper- 
ature. 

• The 3 mm band, which is useful to constrain /3, is less sensi- 
tive to settling than shorter wavelengths even on long baselines (11 
km) and for longer integration times (4 hours). Thus, measuring 
the differential settling between 0.5 and 3 mm would be very time 
consuming. 

• Phase noise should be below about 40° to avoid smearing by 
limited seeing. Although the highest frequencies provide better an- 
gular resolution, this condition favors the 0.8 mm band as the pre- 
ferred frequency to probe the apparent scale height. 

• At the highest inclination (> 85°), the apparent radial depen- 
dency of the surface density is affected by dust settling. However, 
this effect is not a sufficient diagnostic. 

• Other parameters, such as the radial dependency of the dust 
emissivity index, are not substantially altered by settling. 

Our study was performed using a viscosity parameter a ~ 10 -3 . 
Although settling is expected to depend on a, the dependency is 
weak (y/a), and other unknowns, in particular the grain size distri- 
bution but also the surface density, preclude an accurate determina- 
tion of a based on the observation of settling only. 
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APPENDIX A: DETERMINATION OF THE GRAIN 
EMISSIVITY 

Dust grain emissivities can be computed from their dielectric prop- 
erties. However, as grain properties are poorly known in pro to- 
planetary disks, strong assumptions about the grain characteristics 
(shape, composition, porosity, ice layer, ...) have to be made for 
such purpose. 

We follow instead a much simpler approach which takes into 
account the basic asymptotic behaviour of the dust absorption co- 
efficient as a function of wavelength. For wavelengths A much 
smaller than the grain radius a, grains behave as optically thick 
absorbers. Hence, the absorption coefficient per unit mass is wave- 
length independent and 



1000 



v(A, a) — na 2 /m 9 — 3/ (4pda) 



(Al) 



for spherical grains of specific density pa- At the other extreme, for 
A ^> a, the emission coefficient usually falls as: 



«(A,a) oc (l/\f 



(A2) 



where f3 ranges between 1 or 2, depending on the grain composition 
but not on grain sizes. 

In between, for A « 27ra, the absorption coefficient exhibits 
a number of bumps, due to interferences between the refracted and 
diffracted rays. The detailed shape of absorption curve i n this reso- 
nant region depend on grain structure (iNatta et al.ll2004l) . However, 
these detailed shapes will be smeared out when the absorption coef- 
ficient is computed for a size distribution of the grains, so its exact 
knowledge is unimportant provided the overall emissivity curve can 
be reproduced for realistic grain size distributions. 

We thus define the emissivity curve through a small number of 
parameters. For a given grain radius a, the short wavelength regime 
is given by: 



(A3) 



from equation IaTI The long wavelength regime is defined by hi and 
ej, so that for A ^> a, 



K (X) = ki\ ei 

The two regimes intersect at 



A (a) 



k s (a) 



l/ei 



(A4) 



(A5) 



The enhanced emissivity ("bump") is defined at Ai(a) = 
27rha by an enhancement factor f p > 1 compared to the long 
wavelength asymptotic regime 



«(Ai(a), a) = fpki 



Ai(a) 



fpki (27rZi 



(A6) 



The shape around this region is defined by slopes dzeb before (A < 
Ai(a)) and e a after (A > Ai(a)) the bump. The ± sign for e b 
occurs because in this parametrization, «(Ai(a), a) can be smaller 
than the short wavelength asymptotic value k s (a) (see Fig jAll ). The 
emissivity law being a piecewise combination of power laws of A 




Figure Al. Absorption coefficient as a function A for three different size of 
grains(0.1, 1 and 10 mm). These curves correspond to pd = 1.5 g.cm -3 , 
long wavelength parameters exponent e\ = —1.7 and coefficient ki = 0.5 



cm°.g 1 , and bump height f p = 4, exponents e b = 1.5, e a 
and 98 position l\ = 0.65. 



-4.5 




Figure A2. Top panel: Absorption coefficient n at 1 mm as a function of 
CLmax for different exponents p of the size distribution. The solid lines rep- 
resent the results calculated in Ricci et al (2010), using realistic grains made 
of astronomical silicates (10% in volume), carbonaceous materials (20%) 
and water ice. The dashed lines represent our approximated method. Bot- 
tom panel: emissivity exponent (5 computed between 1 and 3 mm for the 
same grain distribution. 



and a, integration over a power law size distribution for the grains 
is straightforward. 

This description with a limited number of parameters cap- 
tures all the required characteristics to adequately represent the 
absorption curves of a given grain size distribution. Figure IA2I 
shows the law used in our s ample models, com pared to the ab- 
sorption coefficients used bv iRicci et al.l J2010bh . The parameters 
are f p = 15.4, e b = 0.68, e a = -3.5, h = 0.65, d = -1.67 
and hi = 0.058. Although differences by 20 % exist, the key fea- 
tures such as the asymptotic values, position width and height of 
the emissivity bump are all well reproduced. 
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